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We study a class of 1+1 quadratically nonlinear water wave equations that combines the linear 
dispersion of the Korteweg-deVries (KdV) equation with the nonlinear/nonlocal dispersion of the 
Camassa-Holm (CH) equation, yet still preserves integrability via the inverse scattering transform 
(1ST) method. This IST-integrable class of equations contains both the KdV equation and the CH 
equation as limiting cases. It arises as the compatibility condition for a second order isospectral 
eigenvalue problem and a first order equation for the evolution of its eigenfunctions. This integrable 
equation is shown to be a shallow water wave equation derived by asymptotic expansion at one 
order higher approximation than KdV. We compare its traveling wave solutions to KdV solitons. 
PACS numbers: 5.45.Yv, ll.lO.Ef, ll.lO.Lm, 47.35. +i 
Keywords: Water-waves, Nonlinear-dynamics, Solitons 



Water wave theory first introduced solitons as solutions 
of unidirectional nonlinear wave equations, obtained via 
asymptotic expansions around simple wave motion of the 
Euler equations for shallow water in a particular Galilean 
frame Later developments identified some of these 
water wave equations as completely integrable Hamilto- 
nian systems solvable by the inverse scattering transform 
(1ST) method, see, e.g., dj. We shall discuss the follow- 
ing 1+1 quadratically nonlinear equation in this class for 
unidirectional water waves with fluid velocity u{x,t), 



mt + cqUx + UTnx + 2mux ^ ~Ju^ 



(1) 



Here m — u — a Uxx is a momentum variable, partial 
derivatives are denoted by subscripts, the constants 
and 7/ Co are squares of length scales, and cq = \/9h 
is the linear wave speed for undisturbed water at rest 
at spatial inflnity, where u and m are taken to vanish. 
(Any constant value u = mq is also a solution.) Equa- 
tion (^ was first derived by using asymptotic expansions 
directly in the Hamiltonian for Euler's equations in the 
shallow water regime and was thereby shown to be bi- 
Hamiltonian and, thus, IST-integrable inJ0]. Before 
classes of integrable equations similar to dfi) were known 
to be derivable from the theory of hereditary symme- 
tries. However, these were not derived physically as 
water wave equations and their solution properties were 
not studied before ||^ . See ||^ for an insightful discussion 
of how the integrable equation (|^) relates to the theory 
of hereditary symmetries. 

The interplay between the local and nonlocal linear 
dispersion in this equation is evident in its phase veloc- 
ity relation, oj/k = (cq — 7 fc'^ ) / ( 1 + a ^ fc^ ) , for waves with 
frequency lo and wave number k linearized around u = 0. 
At low wave numbers, the constant dispersion parameters 
and 7 perform rather similar functions. At high wave 
numbers, however, the parameter properly keeps the 
phase speed of the wave from becoming unbounded. The 



phase speed lies in the band tj/fc € (—7/0;^, co). Longer 
linear waves are the faster provided 7 + cqo^ > 0. 

Equation (Q) is not Galilean invariant. Upon shifting 
the velocity variable by uq and moving into a Galilean 
frame ^ = x — ct with velocity c, so that u(x, t) = 
""(■fj t) + c + Mo, this equation transforms to 



ifit -\- u fh(^ + 2m = — cq — 7 m^^^ 



(2) 



with Co = (cq + 2c + 3uo)i 7 = (7 ~ u^a'^) and appro- 
priately altered boundary conditions at spatial infinity. 
Hence, we must regard equation ([^) as a family of equa- 
tions whose linear dispersion parameters cg, 7 depend 
on the appropriate choice of Galilean frame and bound- 
ary conditions. The parameters cq and 7 may even be 
removed by making such a choice In the following 
we will only use transformations that leave the boundary 
condition m = at spatial infinity invariant. 

This paper reports two main results. First, we identify 
how the dispersion coefficients for the linearized water 
waves appear as parameters in the isospectral problem 
for this IST-integrable equation. We determine how the 
linear dispersion parameters a, cq and 7 in (|l]) affect the 
isospectral content of its soliton solutions and the shape 
of its traveling waves. Second, we rederive equation (|^) 
by using a certain nonlocal form of shallow water wave 
asymptotics that is correct to one order higher than for 
KdV. This new derivation and the analysis we present 
here attaches additional physical meaning to equation 
(|l|) in the context of asymptotics for shallow water wave 
equations. In fact, it even allows us to optimize its dis- 
persion coefficients for maximal accuracy. 

Equation (|l|) restricts to two separately integrable soli- 
ton equations for water waves. When — > this equa- 
tion becomes the Korteweg-deVries (KdV) equation. 
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which for cq = has the famous soliton solution u{x, t) = 
uo sech^((x - ct) yW7/2), c = cq + uq see, e.g., [1). In- 
stead, taking 7 — > in equation (]^) impHes the Camassa- 
Holm (CH) equation, 



ut + cqUx - a Uxxt + 3u Uj; = a 



which for cq = has "peakon" sohton solutions u{x,t) = 
^^-\x-ct\ (discovered and analyzed in 

Equation (|l|) also arises as the Euler-Poincare equation 
for an averaged Lagrangian, , . Its solutions describe 
geodesic motion with respect to the Hi metric of u on 
the Bott-Virasoro Lie group ||^ . The KdV equation arises 
the same way for the L2 metric, see, e.g., 

Equation (|l|) is nonlocal, as is seen by rewriting it as 



-d.r. 



Co 
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and inverting the Hclmholtz operator using its Green's 
function relation (1 — a^d^) e~ ' ^ ' /" = 25{x/a). Drop- 
ping nonlinear terms on the right side and setting 7 = 
recovers the equation studied by Fornberg and Whitham 
1^. Thus, equation contains three simpler water wave 
equations, two of which are known to be integrable. 

Equation ^ is hi-Hamiltonian and hence, isospectral. 
The term bi-Hamiltonian means the equation may be 
written in two compatible Hamiltonian ways, namely as 
mt = -B2{SHi/Sm) = -Bi{SH2/dm) with 

= ^ J u'^ + a ^uj dx, B2 ^ dxm + md^ + 7 9| 

7?2 = ^ J + a'^uu^ ~ ju^ dx , Bi ^ ~ a'^d^ . 

These bi-Hamiltonian restrict properly to those for KdV 
when —>■ 0, and to those for CH when 7^0. Com- 
patibility of Bi and B2 is assured, since d^m + mdx, 
dx and are all mutually compatible Hamiltonian op- 
erators, see, e.g., [0. From this viewpoint, 7 and 
are deformation parameters for the Riemann equation, 
Ut + 3uux = 0. No further deformations of these Hamilto- 
nian operators involving higher order partial derivatives 
would be compatible with B2 [0 . A related approach is 



given in |12 

By the standard Gelfand-Dorfman theory [Q, its bi- 
Hamiltonian property implies that the nonlinear equa- 
tion (|l|) arises as a compatibility condition for two linear 
equations, namely, the isospectral eigenvalue problem. 
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m{x, t) 



and the evolution equation for the eigenfunction "0, 
ipt = -{u + X)il;x + \ux V' ■ 



(4) 



(5) 



Compatibility of these linear equations {il^xxt = 'fptxx) 
and isospectrality {dX/dt = 0) imply equation Con- 
sequently, the nonlinear water wave equation admits 
the 1ST method for the solution of its initial value prob- 
lem, just as the KdV and CH equations do. In fact, 
the isospectral problem for equation (Q) restricts to the 
isospectral problem for KdV (i.e., the Schrodinger equa- 
tion) when ^ and it restricts to the isospectral 
problem for CH discovered in |^ when 7^0. 

The 1ST method for solving the KdV equation was 
introduced in The isospectral problem associated 

with solving the CH equation by 1ST is treated in ||l2|,|l^ . 
The combined eigenvalue problem parameterized by 
a, Cq and 7 is a new variant that retains many of the 
same properties as occur for restrictions to the zero val- 
ues of these parameters. Thus, the 1ST approach based 
on this isospectral eigenvalue problem will illuminate the 
dependence of the solution behavior of equation on 
its dispersion parameters a, cq and 7, in comparison with 
the solution behavior of its two limiting integrable water 
wave equations, KdV and CH. 

Defining a new spectral parameter — J + Xa^ 
yields a spectral problem in the same form as for CH, 



di 



1 



2 / 



7 + c^a^ 
4a2 



which is analyzed and discussed in [n2 17 1. The continu- 
ous spectrum lies in S [0,7-1- cqct). In this form of 
the isospectral equation, the limit to KdV is singular. 

Spectral content and solution behavior. Provided m 
decreases sufficiently rapidly at spatial infinity, equation 
has both continuous and discrete spectra. These spec- 
tral components correspond to the two different types of 
solution behavior available for equation (Q). The con- 
tinuous spectrum of the isospectral eigenvalue problem 
(^) spans the band of allowed linearized phase speeds, 
namely A G (— 7/0;^, co). This continuous spectrum cor- 
responds to radiation (linear waves). The discrete spec- 
trum of lies above this band, with A > cq > 0. The 
discrete spectrum corresponds to the soliton sector of the 
solution space. This is also what is seen in numerical 
computations In the limit that both cq ^ and 

7 — !■ 0, the corresponding isospectral problem for the CH 
equation has purely discrete spectrum representing only 
peakon solutions |^ . 

The derivation from shallow water wave asymptotics is 
similar to that of the KdV equation, except that we keep 
terms of higher order in the small parameters ei = a/h 
and £2 = h'^/P where ei > €2 > and a, h, and I are 
the wave amplitude, the average total water depth, and 
a typical horizontal length scale, respectively. We shall 
sketch the derivation here and give details elsewhere . 
We start with the equations for an inviscid, incompress- 
ible, and irrotational fluid moving under gravity with an 
upper free surface, as, e.g., in |l|]. The velocity potential 
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is expanded at a fixed height zq @ ■ This is preferable to 
the usual depth averaging because it allows a parameter 
freedom in the expansion that may later be optimized. 
The intermediate steps of this calculation all involve the 
chosen height zq. However, the equation for the elevation 
of the surface 77 is the same as in the standard approach 
(see, e.g., ||^). Length is measured in terms of Z, height in 
h and time in I/cq- The elevation rj is scaled with a and 
fluid velocity u is scaled with Coa/h. In nondimensional 
variables one finds 

3 3 1 

= r/f + r/a; + -eiWlx - s^Wvx + ^e2r]xxx + (6) 
Z o D 

/23 5 
+eie2 l^VxVxx + Y^flVxxx 

Transforming to the velocity u gives equation (|^), 
when we expand ry in a series in terms of the fluid veloc- 
ity u = 4>x{x,t; Zq) as 



77 = M + ei (aiw^ + a2Uxd^ ^u) + e2(iuxx 



(7) 



In this transformation we include a crucial nonlocal term 
that was first introduced by Kodama in Usually 
the coefficients of this expansion are constrained by the 
requirement to match the velocity u = 4>x averaged over 
the fluid column. (The derivation in apparently ig- 
nores this constraint.) Considering the fluid velocity at 
height Zq removes this constraint. We introduce another 
optimization parameter by splitting the higher order time 
derivative e2Puxxt obtained by substituting (m) into m) 



PUxxt = -VUxxt-{(i+v) ( Uxxx + '^{iUxUxx + UUxxx) 



which is asymptotically correct to the required order. 
With the choices ai = 2/3 — v, a2 = —5/6 -I- 2v^ and 
/3 — 29/72 — 7z//6 we obtain the equation 

mt+Ux + — (umx + 2mux) + — (1 - &v)Uxxx = , 
2 

where m = u — e2i^Uxx- The parameter u in the splitting 
must be positive in order to obtain correct signs of the 
coefficients. To recover (|^) we reintroduce dimensional 
variables to find = h?v and 7 -I- cqo^ = c^h? /Q. In 
these physical quantities the normalized phase speed is 

kcQ 6 1 -t- vh'^k'^ 

The free parameter v can now be determined by match- 
ing the quartic term in k to the exact result tanh(kh)/kh, 
which gives v — 19/60. The original form of equation (|l|) 
- equation (3) derived in ||] - is recovered by setting 
v = 1/3, which is almost the optimal value. 

Equation (^ determines the fluid velocity u at a fixed 
height Zq. To first order, the velocity u and the eleva- 
tion r] coincide, but the more complicated relation (|^) is 
required to obtain the precise wave elevation profile. 



The traveling wave solution is obtained by the usual 
ansatz u{x,t) = m(s), with s = x — ct, after which (|^) 
can be integrated twice. The solution whose velocity van- 
ishes at spatial infinity is given by 



(c -f 7/a^ — u)a'^{du/ds)'^ = (c — cq — u)u'^ 



(8) 



When the limits of the radiation band coincide at zero, 
the peakon traveling wave equation re-emerges. Oth- 
erwise, the traveling wave solution can be expressed in 
parametric form via a Sundmann transform of the inde- 
pendent variable, (ds/dr)'^ = (c+ j/a^ — u)a^, as 

u{t) — UQsech^ At , A = y/uo/2 , 

sinh At \ 



s{t) = 2a\/ D/uq sinh 



^1-uo/D 



, 1 / D/un — 1 

-2atanh"^ 1 ' ' 



tanh^ At 



-1/2 



where uq — c ~ cq and D — cq + ^ j o? ■ The curvature 
at the maximum is — Mo/2/(7 -I- cqq;^) = — 3u§/(co/i^). 
For a ^ and 7—^0 this smooth family approaches 
the peakon solution. However, this limit cannot be at- 
tained for physical water waves because it implies van- 
ishing mean depth, /i — > 0. In the physical variables the 
solution reads (with rescaled parameter f ) 



u = uosech^^f, A — 



3mo 
2co/i2 



(x - d)A = Vl + o-sinh ^ [^\ + cr sinh Af^ 
— -v/CTtanh^^ ( 1 



1 



a tanh At 



-1/2 



where u = %vuqI cq. The parameter uq/cq plays the same 
role for u as ei = a/h plays for rj. In the KdV limit a —> 
the second equation reduces to x — ct = f. The curva- 
ture at the maximum, — 3mq/(co/i^), is independent of 
u = /h"^ . Increasing v increases the width of the wave 
form without changing the curvature at its maximum. 

Discussion. The water wave equation (^ combines 
the elements of linear dispersion in KdV with the non- 
linear, nonlocal dispersion of CH. This combination is 
IST-integrable, because it retains the bi-Hamiltonian and 
isospectral properties of the KdV and CH equations. The 
dispersion parameters a, co and 7 knit the KdV and CH 
equations into a family of integrable equations and also 
appear together in the eigenvalue problem (^). 

Geometrically, equation (|]) is the Euler-Poincare equa- 
tion for geodesic motion on the Bott-Virasoro group with 
respect to the metric Hi as shown in |^ . One may choose 
an arbitrary metric for the Euler-Poincare equation ex- 
pressing geodesic motion on either the diffeomorphism 
group or the Bott-Virasoro group, and again seek condi- 
tions for integrability. Integrable cases of this geodesic 
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motion for norms other than L2 and Hi are not yet 
known, although the elastic scattering properties of the 
confined pulses in some of these other cases are discussed 
in pl[ |. In particular, the Galilean invariant integrable 
limiting equation, 

'^xxt ~^ '^'^x^xx ~t~ '^^xxx — n^x ^'^xxx : (9) 

a 

is obtained by taking — s- cx) in (|l|) while the ratios 
co/a^ and 7/0;^ remain constant. For cq and 7 absent, 
this equation has soliton solutions with compact support 
(compactons) that were studied numerically in an- 
alytically in |2^] on the real line, and analytically in |2^] 
for the periodic case. 

The dispersion parameters cq and 7 in (|l|) represent 
combined velocity shifts and Galilean boosts. These 
transformations are deformations of the CH equation - 
they are not its symmetries. While these dispersive defor- 
mations preserve the isospectral property of the CH equa- 
tion, they alter its isospectral content in two ways. First, 
they add to the discrete CH peakon spectrum a band 
of continuous spectrum near the origin. This isospectral 
band corresponds to linear radiation and spans the al- 
lowed range of linear phase velocity. Second, either cq 
or 7 nonvanishing breaks the reflection-reversal symme- 
try of the CH equation for Cq = that allows coexistence 
of its generalized-function peakon and antipeakon soliton 
solutions traveling in opposite directions. Thus, adding 
linear dispersion breaks this reflection-reversal symme- 
try and removes the antipeakon solutions, while it also 
rounds the peaks of the peakon solutions, thereby regu- 
larizing them into ordinary solitons, plus linear radiation. 

This letter reconfirmed that is a genuine shallow 
water equation by rederiving it via asymptotic expansion 
of the Euler equations. This derivation used a nonlocal 
transformation due to Kodama as well as a splitting 
of the time derivative at the highest order that introduced 
an optimization parameter, ly — /h? . We showed that 
the equation originally discovered in |^ via asymptotic 
expansion of the Eulcr fluid Hamiltonian is contained in 
the family of equations derived here for the case v = 1/3. 
Moreover, the value 1/ = 1/3 obtained in the Hamiltonian 
derivation in differs only very little from the optimal 
value V = 19/60 yielding the most accurate linear dis- 
persion relation among those in the family derived here. 
We also showed that varying v alters the width of the 
traveling wave proflle for equation (|l|) without changing 
the curvature at its maximum. 
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